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Are  obtained  equations,  which  describe  laminar  flame  in 
heterogeneous  hydrodynamic  field  if  curvature  of  flame  front  can  be 
disregarded  (action  of  large-scale  turbulence).  The  basic  properties 
of  this  flame  are  clarified.  The  possible  mechanism  of  extinguishing 
of  flame  is  obtained.  The  mechanism  of  turbulent  combustion  is 
considered.  t ^  -A  I  ^  ° 


(*?>  y  /«  * 

J',  \ 


i 

i 

1.  Heterogeneity  of  hydrodynamic  field  has  fundamental 
importance  for  mechanism  of  turbulent  combustion;  we  will  consider 
turbulence  as  probability  field  of  velocity  gradients.  Let  us 
designate  through  i  the  size  of  the  region,  in  which  the  velocity 
gradients  for  a  certain  period  of  time  (lifetime  of  this  region) 
weakly  depend  on  the  position  of  point  within  the  region,  although 
they  can  arbitrarily  depend  on  time.  This  size  is  identical  to  the 
scale  of  turbulence.  Let  us  introduce  also  the  significant  dimension 
5  of  "temperature"  field,  after  defining  it  as  distance  in  the 
direction  of  the  gradient  of  the  temperature  (or  concentration),  on 
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which  the  temperature  (concentration)  substantially  varies.  In  our 
case  this  will  be  the  thickness  of  laminar  flame  in  the  homogeneous 
mixture.  The  character  of  phenomenon,  obviously,  strongly  depends  on 
the  relationship  between  l  and  6.  If  l» 6  (large-scale  turbulence), 
flame  front  weakly  is  bent  (r>>6,  where  r  -  radius  of  curvature  of 
flame  front),  but  the  area  of  its  surface  continuously  is  changed: 
flame  is  stretched  or,  on  the  contrary,  is  contracted.  Is 
investigated  below  precisely  this  case  (/>> 6). 

Let  us  consider  field  of  velocities  of  region  l  without 
combustion,  which  we  will  call  external  hydrodynamic  field.  The  field 
of  pure  shear  (Fig.  la)  is  the  simplest  example,  when  \e  correlation 
between  the  velocity  components  is  absent.  Another  example  is  "flow 
at  critical  point”  (Fig.  lb),  that  appears,  if,  for  example,  in 
regions  A  and  B  rarefactions  arose.  The  correlation  between  the 
velocity  components  here  is  exhibited  in  the  pure  form.  The  velocity 
field  of  arbitrary  region  /  is  a  superposition  of  these  two  limiting 
cases.  For  example,  in  the  plane  case,  taking  into  account  the 
equation  of  continuity,  velocity  field  takes  the  form 

u"  =  k  (t)r  b,  (/)  y.  v°  —  —  k(t)  y  MO-r 

It  is  evident  that  to  "flow  at  critical  point"  u=kx,  v=-ky  are 
superimposed  arbitrary  shifts  u=k,y,  v=k2x.  The  three-dimensional 
case  in  principle  does  not  differ  from  the  plane. 
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Fig.  la,  b. 


Page  50. 

Assumptions,  essential  for  obtaining  equations,  will  be  usual 
assumptions  of  theory  of  slow  combustion:  1)  pressure  constancy 
(equation  of  motion  is  not  examined),  2)  in  equation  of  energy  is 
considered  only  thermal  energy  of  gas. 


Furthermore,  we  will  consider  only  primary  processes  of  transfer 
and  consider  that  heat  release  rate  is  determined  by  concentration  of 
one  component  of  gas  mixture  and  by  temperature,  and  Lewis  number  is 
equal  to  1.  Then,  as  is  known,  the  equation  of  diffusion  it  is 
possible  to  exclude. 


In  conventional  designations  in  plane  case  we  will  have 
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Furthermore,  is  assigned  applied  field  u°,  v°.  Let  the  applied 
field  -  pure  shear.  Without  the  limitation  of  generality  it  is 
possible  to  consider  u°=a(t)y,  v°=0;  the  orientation  of  the  initial 
distribution  of  temperature  (being  one-dimensional)  relative  to 
applied  field  is  assigned  by  angle  vP,  where  <p<p  -  flame  front  (Fig.  2). 

True  hydrodynamic  field  will  differ  from  external,  since  during 
combustion,  as  a  result  of  thermal  expansion  of  gases,  motion,  normal 
to  flame  front,  appears.  As  a  whole  the  phenomenon  is  such:  flame 
front  is  turned  by  applied  field  in  plane  xy  with  a  change  in  the  area 
of  its  surface,  simultaneously  it  moves  relative  to  gas,  creating 
additional  velocity  field  due  to  the  thermal  expansion. 

We  seek  solution  of  equations  (1.1)  in  the  form 

f 

T  =  T  (t,  C).  u  =  ay  -  t  U.  ;>[l  -  (iff*  —  \adtX  1 

0 

f 

r  =  v  (t.  ;)[tg*  —  W<][1  +  ( tg t  -W)] 

0  0 

where  f  , 

;  =  [*  -  »(tgi|:  -$«*)][!  T  (tpt  -$«&)*]'’  ’ 

0  « 


For  determining  T(t,  $)  and  V(t,  $)  we  will  obtain  equations 


<*T 

f'flv, 


,  ,  ,  dT  1  »  9T\  r  T 

f  -  1  r  s;)~F  {T) 


*  '  <  t  £ 


(1.2) 
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Here 


<p  (/)  =  a  [tg  —  \arf/|  1  +(lgif 


\a*)T 

0 


rdinate  $  is  always  normal  to  the  surface  of  flame. 
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Fig.  2. 


Page  51. 


Velocity  V  is  connected  with  expansion  of  gas.  If  combustion 
products  and  combustible  mixture  occupy  half-spaces,  it  is  natural  to 
consider  V=0  for  the  fresh  mixture.  Initial  and  boundary  conditions 
for  equations  (1.2)  in  this  case  will  be 


7  -  T  (l)  when  t  =  o,  v  _  o  with 


If  at  initial  moment  of  time  combustion  products  occupy  the  band 
of  finite  width  with  the  symmetrical  temperature  distribution,  the 
simplest  initial  and  boundary  conditions  will  be 


T=rC)  u>K«^  f  =  o, 


^  =  v  =  o  w  ;  =  o 


Let  now  u°=-k(t)x,  v°=k(t)y.  We  seek  the  solution  of  equations 


(1.1)  in  the  form 


T  —  T  (t,  S) 


i  * 

u  =  -  Ax  4-  V  (t,  o  ag*  exp(  2  [kdt'j  [l  4-  ctg*^  exp  (4^  k  rf/)  j“  * 

0  0 

f  t 

r  =  ky  +  V  {t,  5)[i  +  ctgH  exp(4^/r*jp 
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where 

i  t  t 

;  =  [y  +  i  ctgifexp  +  ctg^exp^Ad^J”'’ 

0  0 

For  determining  T(t,  $)  and  V(t,  $ )  we  will  obtain  (1.2),  where 

i  i  _j 

q  (/)  =  A- [l  -  ctgH  exp(4^*dt)][l  4-  ctg* exp(4^**)] 

L  ft  0 

In  general  case  we  will  also  obtain  equations  (1.2).  The 
specific  character  of  each  case  is  included  in  function  <p(t). 

It  is  not  difficult  to  check  that  «p(t)  is  relative  rate  of  change 
of  surface  area  of  flame,  taken  with  opposite  sign. 

Equations  (1.2)  it  is  possible  to  also  obtain,  having  assumed 
that  flame  remains  plane  and  is  assigned  function  <p(t).  Let  £ ,  r? ,  $  - 
coordinate  system,  oriented  so  that  the  axis  $  is  always  normal  to  the 
surface  of  flame,  and  let  vt,  -  corresponding  velocities. 


Since  surface  of  flame  is  plane, 


Consequently , 


(1-3) 
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Let  us  consider  small  element  of  surface  d£  di?  in  arbitrary  plane 
$=const.  Relative  rate  of  change  in  its  area 


1  d(dldr])_^jt 

-?(')-  dl  d\  ^  dr\ 


(1.4) 


Page  52. 

After  representing  r-  in  the  form  l;  ((,  ()  =  n  (/)  ;  V  (t,  l)  and  taking 
into  account  (1.4),  we  will  obtain  from  (1.3)  equation  (1.2). 

Let  us  introduce  dimensionless  quantities 

t  =  tu0/60,  f>  =  t/60.  «o  =  r/u0,  o  =  (r  -  r0) / (r,  -  t0) 

p'  —  P  /  p0.  *■'  —  >•  i  K<  /  =  F  '  F0,  £  /  a0  (6a  =  >.)  /  \ 

Po  =  P  (T'o).  =  x  (T0),  F„  =  pflU.f  p  (Tl  -  T0)  /  d0 

Here  u0  -  velocity  of  propagation  of  normal  flame;  60  -  distance, 
at  which  in  normal  flame  8  is  changed  approximately  e  times  (exactly  e 
times  in  solution  of  Michealson) ;  T0,  T,  -  initial  and  final 
(adiabatic)  temperatures  of  gas;  F0  -  average  heat  release  rate  in 
normal  flame. 

Parameter  e  is  ratio  of  characteristic  time  of  smooth  burning 
t0  =  60/u0  to  characteristic  time  of  action  tl  =  l/(p  on  flame  of  external 
hydrodynamic  field;  ti  -  time,  for  which  with  <p=const  surface  area  of 
flame  is  changed  e  times,  or,  which  is  the  same,  time,  within  which 
distance  between  surfaces  of  constant  temperature  would  be  changed  e 
times  in  the  absence  of  chemical  reaction  and  thermal  conductivity. 
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In  adopted  designations  system  (1.2)  takes  form 


p 


3 ^ 
d3 


d_ 

as 1 


(>■1;+ /<♦> 


S+'<T)Pg  + 


a  (po.  >  n 
03 


(1.5) 


Here  and  throughout  primes  with  p  and  A  are  omitted. 


New,  in  comparison  with  equations  of  one-dimensional  flame,  terms 
in  (1.5)  are  of  order  e(r).  Let  us  do  numerical  estimations.  Let, 
for  example,  the  velocity  of  the  turbulent  pulsations  v'~10  m/s  with 
the  scale  of  turbulence  /~1  cm;  then  t^lO"3  s.  For  the  majority  of 
the  fuel-air  mixtures,  close  to  stoichiometric,  to~10'3  s  and  rapidly 
increases  to  10‘ 1  s  during  impoverishment  or  enrichment  of  mixture 
approximately  1.5  times.  Respectively,  we  will  obtain  e~l-10, 
depending  on  the  composition  of  mixture. 


Parameter  e  is  in  close  agreement  with  parameter  7=v'80/u0/, 
proposed  [1]  for  evaluation  of  relative  role  of  "surface"  and 
"volumetric"  combustion.  The  difference  between  them  is  that  in  7 
value  v'/l  is  the  averaged  velocity  gradient,  and  value  <p  in  c  - 
instantaneous  and  local.  Therefore  e  can  be  both  positive  and 
negative  value.  It  is  evident  from  (1.4)  that  the  gradients  of 
different  velocity  components  can  make  different  contribution  to  a  : 

change  in  the  surface  area  of  flame  -  flame  can  be  stretched  or 
contracted  in  all  directions,  but  it  can  be  contracted  in  one 

direction  and  be  stretched  in  another.  • 

I 

1 

i 

1 

I 
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Let  us  write  equations  (1.5)  in  coordinate  system,  connected  with 
flame  front,  i.e.,  in  system  of  coordinates,  whose  origin  is  connected 
with  certain  temperature  de.  The  coordinate  /30  of  point  with  the 
temperature  d0  is  determined  by  the  equation 


dl°  =  [e  (t)  fi0  —  —  |m  ft)  -f  to]  -  £  (t)  m  (t) 

dx 

I 

Consequently, 


T  t  •  T 

■30  (t)  =  exp  {  e  (t)  dx  jpoo  —  \  m  (r)  exp  \  t  (x)  dx  j]  (3ud  =  3„(0)) 


Page  53. 

In  variables  a=/3-/3„,  t  system  (1.5)  will  take  form 


)V> 

[""t 


•  ta  +®)Sl==-<sr(x  ~L )~  1  ({'] 


t»T 


(m  +  €a) 


d  (puO  _  q 
eta 


(1.6) 


Function  m(r)  is  determined  during  solution  of  specific  problem. 
The  initial  value  m(r)  must  be  assigned. 


Equations  (1.6)  can  be  also  obtained  from  equations  (1.5),  by 
considering  that  latter  are  obtained  from  (1.3),  where  £,  r?,  $  - 
coordinate  system,  connected  with  flame  front.  It  is  necessary  only 

to  represent  cj  in  the  form 

u,  (3.  x)  =-  m(T)  4- <0,(3,  x).  m(x)  --  “>(-<*•  x> 
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2.  Before  constructing  solution  of  problem,  let  us  give 
qualitative  description  of  picture  of  phenomenon  in  question. 

In  uniform  hydrodynamic  field  (e=0)  solution  of  system  (1.6) 
under  some  assumptions  relative  to  function  of  heat  release  f(fl)  and 
corresponding  initial  and  boundary  conditions  approaches 
time- independent  solution  and  which  describes  temperature  distribution 
in  normal  flame.  In  this  case  m(T)-*l.  The  stationary  form  of  normal 
flame  is  a  result  of  the  equilibrium  of  opposite  "forces":  the 
smoothing  out  action  of  thermal  conductivity  is  balanced  by  the 
opposite  action  of  chemical  reaction. 

With  e *0  appears  new  factor  -  convective  heat  transfer  ea 
which  either  smooths  out  (e>0),  or  it  increases  (e<0)  slope  of 
temperature  distribution. 

If  e-*const<0,  solution  of  system  (1.6)  must  approach  certain, 
steeper  than  normal  flame,  stationary  temperature  distribution. 

Different  proves  to  be  the  situation  with  e>0.  In  the  model  of 
gaseous  medium  accepted  the  thermal  disturbances  are  propagated  with 
infinite  velocity.  In  practice  this  velocity  is  also  very  great 
(speed  of  sound).  Therefore  the  thermal  effect  of  hot  region 
instantly  penetrates  to  large  distance  into  the  cold  region.  At  the 
same  time  at  the  sufficiently  large  distance  from  the  origin  of 
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coordinates  to  the  side  of  cold  gas  (i.e.  with  the  sufficiently  large 
negative  a)  the  convection  heat  current  acquires  the  same  direction  as 
the  heat  flow,  caused  by  thermal  conductivity.  During  the  slow 
combustion  such  heat  transfer  into  the  cold  region  is  compensated  in 
no  way;  conditions  at  the  large  distance  from  the  flame  front  always 
vary  and  entire  process  becomes  unsteady. 

For  quantity  q  of  heat  seperating  per  unit  of  surface  of  flame 
per  unit  time  (proportional  to  quantity  of  burning  substance),  taking 
into  account  that  with  a -»±oe,  relationship 

<1  -  \  /da  (  p  ^  -  (  p  (rn  -  ex  u>)  ^  dx 

-  m  — oc  -cc 

In  normal  flame  q-p (m+w) =const .  In  remaining  cases  (e*0  or  e=0, 
but  phenomenon  is  unsteady)  q  is  not  connected  unambiguously  with 
p(m+u)  and,  generally  speaking,  it  is  variable. 

In  time  intervals,  when  applied  field  vanishes,  flame  is 
one-dimensional,  but  unsteady. 

Structure  of  this  flame  has  much  in  common  with  flame  structure 
with  e*0;  therefore  let  us  consider  it  in  more  detail. 

Page  54. 


In  one-dimensional  unsteady  flame  m+a>,  being  velocity  of  motion 
of  temperature  d0  relative  to  gas  (m  -  relative  to  fresh  mixture,  m+u> 
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-  relative  to  particles  of  gas  vith  other  temperatures),  can  take  any 
values,  in  particular  negative.  If  m+u<0,  gas  flows  into  this  plane 
a-const  not  from  the  side  of  fresh  mixture,  but  from  the  side  of 
combustion  products.  The  flow  of  molecules  (as  heat  flow)  is  composed 
of  the  flows  of  convective  and  diffusion.  In  the  normal  flame 
convective  and  diffusion  flows  of  fuel  are  combined.  In  the  unsteady 
flame  they  can  be  subtracted.  Selecting  conditionally  the  zone  of 
reaction,  it  is  possible  to  encounter  the  situation,  when  gas  flows 
into  this  zone  from  the  side  of  combustion  products,  and  fuel  is 
supplied  towards  the  convection  current  by  diffusion.  Let,  for 
example,  initial  function  for  equation  (1.5)  be 

d(3,0)  =  0  when  3<o,  ^ (3.0)  =  1  with  3>°  (2-*) 

before  will  be  established  the  normal  flame,  in  which  any 
temperature  moves  towards  the  fresh  mixture,  for  a  certain  period  of 
time  the  points  with  the  temperature,  close  to  d=l,  will  move  towards 
the  combustion  products.  Consequently,  at  least  in  the  part  of  the 
zone  of  reaction,  the  convection  current  of  fuel  (in  the  system  of 
coordinates  a)  will  be  directed  against  the  diffusion.  The 
establishment  of  normal  flame  in  this  case  occurs  by  means  of  the 
leveling  action  of  thermal  conductivity  and  diffusion. 

If  function  d(/3,  0),  on  the  contrary,  is  very  flat  in  comparison 
with  case  of  normal  flame,  main  factor  of  establishment  of  normal 
flame  will  be  local  increase  in  temperature  due  to  chemical  reaction. 
In  this  case  the  velocity  of  motion  relative  to  gas  of  points  with 
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temperatures,  close  to  d=l,  at  the  first  moments  of  time  will  be 
considerably  more  than  this  velocity  in  the  normal  flame,  having  the 
same  direction. 


It  follows  from  equations  (1.6)  that  in  heterogeneous 
hydrodynamic  field  the  concept  of  velocity  of  propagation  of  flame 
front,  strictly  speaking,  does  not  make  sense.  With  any  the 
velocity  of  the  motion  of  gas  relative  to  the  fixed  temperature 
varies  in  the  range  (-»,  +1*)  ,  depending  on  the  temperature  of  the 
particles  of  the  gas.  Certainly,  with  small  |  e  |  this  concept  can 
make  approximate  sense.  Hydrodynamic  field  in  the  system  of 
coordinates  of  flame  appears  on  "flow  at  critical  point"  mentioned 
above  with  the  plane  of  symmetry,  parallel  to  flame  front.  With  small 
|  c  |  this  plane  is  located  far  from  the  region  with  the  high 
gradients  of  temperature  (region  of  basic  part  of  the  change  in  the 
temperature  from  zero  to  one);  the  field  of  the  mass  flow  rate  in  the 
latter  is  close  to  uniform  and  flame  differs  little  from  the 
one-dimensional.  With  the  sufficiently  large  |  e  |  the  plane  of 
symmetry  is  moved  into  the  region  with  the  high  temperature  gradients 
and  the  concept  of  the  velocity  of  propagation  of  flame  front  does  not 
have  even  approximate  sense. 

3.  Integration  of  first  equation  (1.6)  with  respect  to  t  gives 


$  (a,  t)  =  (a,  0)  -  \  (m  -  ea  +  <■>) 


If  flame  was  normal  at  initial  moment,  for  small  r  with  e=const 
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we  will  obtain 

0  (a,  t)  =  6  (a,  0)  —  eat  -f  0  (r*) 


Page  55. 

In  normal  flame  dd/da20;  consequently,  temperature  distribution 
begins  to  be  smoothed  out  with  values  e>0  and  to  become  steeper  with 
e<0 . 


We  investigate  behavior  of  flame  with  large  |e|.  In  this  case 
the  temperature  distribution  everywhere,  with  exception  of  the  narrow 
layer  of  thickness  h~l/VT*  I  in  the  vicinity  of  point  a=0,  is 
described  by  the  equation 


i  cfd  .  d&  n 
—  _  -f  a  —  =  0 
•  dx  da 


(3.1) 


which  is  obtained  from 

I 


(1.6)  by  the  rejection  of  small  terms. 


This  equation  possesses  system  of  characteristics  (x=aexp(  er ) , 
along  which  G^const. 

Let  e<0.  Then  it  follows  from  (3.1)  that  for  the  time  r~l/|e  | 
almost  all  the  change  in  the  temperature  will  be  concentrated  in  layer 
h.  With  |  e  |  ■*<*>  we  have  h->0,  and  the  temperature  profile  approaches 
disruptive  1 . 


FOOTNOTE  1 .  If  we  take  into  account  thermal  conductivity,  for  the 
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temperature  distribution  we  will  obtain 

0  =  y  +  \  erf  ( (—  oc  <  a  <  +  oc ) 

passing  with  e-*-®  in  (2.1),  if  we  in  (2.1)  by  /3  imply  a.  ENDFOOTNOTE. 

With  c>0  (3.1)  gives  "deliquescence"  of  flame  with  the  same 
characteristic  time  T~l/e;  distance  between  points  with  any  two 
temperatures  grows  proportionally  to  exp  (ft).  respectively  everywhere 
decreases  gradient  of  temperature. 
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Hence 


where  6#Sdlf  d2,  d2,  d4Sl.  In  obtaining  (3.4)  was  used  the 
integration  by  parts  and  was  used  condition  36/9a=0  with  6=1. 

Function  a>  (gas  velocity,  connected  with  thermal  expansion)  must 
have  limited  derivative  9cj/9a. 

Page  56. 

Completing  in  (3.4)  passage  to  the  limit  with  6„-»l,  we  will 
obtain  (3.3).  Let  us  note  that  in  this  case 

q  —  const  =  X  (1)  j 

Thus,  using  function  (3.2),  we  will  have  6=1  with  a£0, 

96/9a=const  with  a=-0.  In  the  field  [-<*,  0)  the  function  6  is 
described  by  equations  (1.6)  with  f=0. 

We  will  seek  steady-state  solution  of  these  equations  with 
e=-c=const<0 .  For  obtaining  the  qualitative  picture  we  will  be 
bounded  to  approach  p=X=l.  Taking  into  account  that  96/9oc=l  with  a=0, 
let  us  find 

a 

ft  =  \  e-*1'  **"“  da 


-oc 
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Value  m  is  determined  from  condition 

0 

^  da  =  1 

-CO 

In  particular  h,  6 

m  —  \  npn  c  =  0,  m  =  0  npn  c  =  nl  2 
m - \^2c  In  \  "c  npn  e  -  oc 

Key :  ( 1 ) .  with . 

Solutions  for  different  c  are  given  in  Fig.  3. 

With  c>jt/2  we  have  m<0,  i.e.,  gas  flows  into  zone  of  reaction 
(point  a=0)  from  side  of  combustion  products.  Curve  Ma)  at  point 
a=-0  has  positive  curvature  with  m>0  and  negative  with  m<0.  On  the 
section  of  convexity  of  this  curve  the  gas  velocity  m=ca  is  negative 
and  diffusion  flow  of  fuel  is  directed  against  the  convective. 

Sequence  of  curves  d(a)  with  change  of  c  from  ®  to  0  is 
qualitatively  identical  to  sequence  of  curves  during  establishment  of 
normal  flame  from  disruptive  profile  of  temperature  (2.1). 

With  c>0  steady-state  solution  of  equations  (1.6),  which 
satisfies  conditions  0{-®)=O,  d(+®)=l,  does  not  exist.  The  physical 
causes  for  this  were  explained  in  p.  2. 


However,  with  e<<l  it  is  possible  to  obtain  approximate 
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steady-state  solution,  after  introducing  fictitious  heat  of  sink,  with 
d*0.  If  the  power  of  the  heat  of  sink  is  small  in  comparison  with  the 
power  of  heat  source,  it  is  possible  to  consider  that  the  temperature 
distribution  will  evolve  slowly  and  for  a  certain  period  of  time  will 
differ  little  from  the  solution  with  the  sink. 

Using  (3.2),  with  p=X=l  in  stationary  case  we  will  obtain 

Heat  sink  q0  will  be  minimum,  if  it  is  located  at  point  a„=-m/e 
(q0=exp(-m2/2e ) ) .  Boundary  conditions  for  (3.5) 

ft  ^  i  when  3  ",  ft  =  0  with  >  = 

Consequently , 

2 

ft  ^  i  rf(l  (- 

-  rr.  t 

Condition  for  determining  m  will  be 

o 

\  *  da  =  1  (,'■!.()) 

~m  t 

Page  57. 

During  solution  of  (3.6)  relative  to  m  appears  ambiguity,  easily, 
however,  removable  by  condition  m=l  with  e=0. 
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We  give  results  of  calculations  of  m  and  aB  for  some  values  of  e 

e  =  0  0.097  0. 18  0.26  0.40  0.516  0.57  0.58 

m  1.0  1 .in  1.20  1.26  1,26  1.22  1.07  1.0 

—  Oo=  oo  11.3  6.67  4.85  3.16  2.36  1.88  1.72 

Curves  d(a)  are  given  in  Fig.  3.  Close  to  reality  will  be  the 
solutions,  apparently,  only  to  e*=0.2,  when  the  value  of  sink  is  still 
small  (q0=exp  (-2)  with  e=0.18).  With  e>0.3  the  effect  of  sink  is  so 
considerable  that  m  begins  to  decrease.  The  solution  does  not  exist 
with  f>0.58  (with  e=0.58  we  have  qc*l/e).  Thus,  with  e~l  phenomenon 
becomes  essentially  unsteady. 

4.  Use  of  function  (3.2)  does  not  make  it  possible  to  find 
change  of  q.  For  solving  this  question  it  is  necessary  to  use  zone  of 
reaction  of  finite  width.  Nevertheless,  the  explained  general 
properties  of  phenomenon  make  it  possible  to  draw  qualitative 
conclusions  without  the  exact  solution  of  problem.  Actually 

■too  1 

-oo  0 

With  e <0  in  region  of  basic  part  of  change  of  temperature  (and, 
consequently,  in  basic  part  of  zone  of  reaction)  da/d-d  is  reduced  in 
comparison  with  normal  flame,  with  e>0  -  it  increases.  In  other 
words,  is  narrowed  down  or  is  widened  the  zone  of  reaction; 
respectively,  q  is  reduced  or  increased. 

It  is  obvious  that  effect  of  e  on  q  will  be  greater,  the  wider 
the  zone  of  reaction  in  normal  flame. 
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If  we  consider  quantity  Q  of  heat,  which  is  isolated  per  unit  of 
primary  surface  of  flame,  then  it  is  possible  to  show  that  with  time  Q 
begins  to  increase  with  e<0  and  to  be  reduced  with  e>0,  if  at  initial 
moment  the  flame  was  normal.  With  e<0  in  the  stationary  case 
Q~qexp(cr).  In  the  turbulent  flow  the  stretching  of  flame  (e<0) 
statistically  predominates  above  contraction  (e>0).  This  directly 
follows  from  the  form  of  the  function  <p(t),  but  it  can  be  considered 
also  as  experimental  fact.  Result  -  increase  in  rate  of  combustion  in 
the  presence  of  turbulence. 

In  flame  with  chain  mechanism  of  chemical  reaction  the  change  of 
thickness  of  flame  will  influence  diffusion  of  active  centers  and, 
apparently,  q  will  be  changed  more  slowly  than  with  thermal  mechanism 
of  chemical  reaction. 

Is  possible  also  change  in  reaction  mechanism:  decrease  of 
thickness  of  flame  will  shift  mechanism  in  direction  of  chain, 
increase  -  in  direction  of  thermal. 

5.  If  f({>)=0  with  d<£o>0,  f(d)2:0  with  $,,£$<1,  then  with 
sufficiently  large  c  there  is  no  different  from  zero  steady-state 
solution  of  equations  (1.5),  which  satisfies  boundary  conditions 


d-*0  when  3^±oc, 


-  =  w  =  0  with  9  =  o 
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Page  58. 

I.e.  with  sufficiently  strong  and  prolonged  elongation  finite 
volumes  of  combustion  products,  surrounded  by  fresh  mixture,  can 
vanish.  The  process  of  such  extinguishing  of  flame  is  very  simple. 

Let  at  initial  moment  of  time  the  combustion  products  occupy  in  the 
plane  $  the  band  of  finite  width  d>>60.  With  sufficiently  large  c  the 
contraction  of  band  under  the  action  of  external  hydrodynamic  field 
cannot  be  compensated  by  the  process  of  combustion  and  the  bandwidth 
begins  to  be  reduced.  If  in  this  case  is  reached  d~60,  the  effect  of 
cold  regions  is  closed  and  occurs  the  temperature  decrease, 
accompanied  by  end  of  burning. 

6.  Investigation  conducted  above  makes  it  possible  to  make  some 
conclusions  about  mechanism  of  turbulent  combustion  (in  the  case  of 
large-scale  turbulence). 

Important  role  plays  parameter  7=v'60/u0/.  Order  of  magnitude  7, 
obviously,  coincides  with  the  probable  value  of  the  absolute  value  of 
parameter  e=<p60/u0. 
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With  7<<1  combustion  occurs  in  laminar  fronts.  Perceptible 
deviations  of  flame  from  the  normal  and  the  phenomena,  connected  with 
extinguishing  of  flame  (resorption  of  finite  volumes  of  combustion 
products,  discontinuities  of  flame,  the  throwing  of  combustion 
products  into  the  combustible  mixture),  are  unlikely. 

With  7~1  structure  of  laminar  flames,  as  a  rule,  differs  from 
normal:  concept  of  velocity  of  propagation  of  flame  front  becomes 
meaningless,  temperature  profile  and  quantity  (specific)  of  burning 
substance  are  variable.  The  probability  of  extinguishing  of  flame 
increases;  but,  taking  into  account  that  for  the  possibility  of 
extinguishing  of  flame  |  e  |  must  be  sufficiently  large  (c>ir/2  in  the 
idealized  diagram  of  p.  3)  and  is  necessary  some  time,  the  greater  it 
is,  the  less  |  e | ,  this  probability  should  be  recognized  still  small. 
Thus,  with  7~1  combustion  occurs  predominantly  in  the  "frontal"  mode; 
the  contribution  of  the  "volumetric"  combustion,  connected  with 
extinguishing  of  flame,  is  small. 

With  7>>1  structure  of  laminar  flames  strongly  differs  from 
normal,  probability  of  extinguishing  of  flame  becomes  considerable. 
Combustion  occurs  both  in  laminar  fronts,  which  constantly  burst  and 
extinguish,  leaving  "warm"  regions,  and  with  the  help  of  volume 
reactions,  which  go  in  "warm"  regions  and  constantly  are  accompanied 
by  ignition  with  formation  of  laminar  flame. 
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